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Abstract 

The rules associated with prepositional logic programs and the sta- 
ble model semantics are not expressive enough to let one write concise 
programs. This problem is alleviated by introducing some new types 
of prepositional rules. Together with a decision procedure that has 
been used as a base for an efficient implementation, the new rules sup- 
plant the standard ones in practical applications of the stable model 
semantics. 



1 Introduction 

Logic programming with the stable model semantics has emerged as a viable 
method for solving constraint satisfaction problems [4, 5]. The state-of- 
the-art system smodels [6] can often handle non-stratified programs with 
tens of thousands of rules. However, prepositional logic programs can not 
compactly encode several types of constraints. For example, expressing the 
subsets of size k of an n-sized set as stable models requires on the order of 
nk rules. In order to remedy this problem, we improve upon the techniques 
of smodels, by extending the semantics with some new types of propositional 
rules: 

• choice rules for encoding subsets of a set, 

• constraint rules for enforcing cardinality limits on the subsets, and 
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• weight rules for writing inequalities over weighted linear sums. 

The extended semantics is not based on subset-minimal models as is the 
case for disjunctive logic programs. For instance, the choice rule is more of 
a generalization of the disjunctive rule of the possible model semantics [7]. 

A system that computes the stable models of programs containing the 
new rules has been implemented [9], and it has successfully been applied 
to deadlock and reachability problems in a class of Petri nets [3]. Other 
problem domains, such as planning and configuration, will benefit by the 
improved rules as well. The system is based on smodels 1.10 from which it 
evolved. 

The new rules and the stable model semantics are introduced in Sec- 
tion 2. A decision procedure for the extended syntax is presented in Sec- 
tion 3, and some important implementation details arc described in Sec- 
tion 4. Experimental results are found in Section 5. Readers not familiar 
with monotonic functions should consult the appendix. 

2 The Stable Model Semantics 

Let Atoms be a set of primitive propositions, or atoms, and consider logic 
programs consisting of rules of the form 

h <— ai, . . . , Qji, not bi, . . . , not bm, 

where the head h and the atoms oi, . . . , a„, 6i, . . . , h„i in the body arc mem- 
bers of Atoms. Call the expression not b a not-atom — atoms and not-atoms 
are referred to as literals. 

The stable model semantics for a logic program P is defined as follows [2] . 
The reduct of P with respect to the set of atoms A is obtained by 

1. deleting each rule in P that has a not-atom not x in its body such 
that X E A, and by 

2. deleting all not-atoms in the remaining rules. 

Definition 1. A set of atoms 5 is a stable model of P if and only if S is 
the deductive closure of P^ when the rules in P^ are seen as inference rules. 

In order to facilitate the definition of more general forms of rules, we 
introduce an equivalent characterization of the stable model semantics. 
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Proposition 1. We say that fp : 2^*°™* ^ 2^*°™* is a closure if 

fp{S) = {/i I /i <— ai, . . . , an, not bi, . . . , not bm G P, 

oi, • • • , On e fp{S), bi,...,bm^ S}. 

Let 

gp{S) = P|{/p(5) I fp : 2^*°™^ ^ 2^*°"^ is a cZosure}. 

Then, S is a stable model of the program P if and only if 

S = gp{S). 

Proof. Note that the deductive closure of the reduct P^ is a closure, and 
note that for every fp that is a closure, the deductive closure of P^ is a 
subset of fp{S). □ 

A stable model is therefore a model that follows from itself by means of 
the smallest possible closure. In other words, a stable model is a supported 
model, and this is the essence of the semantics. 

Definition 2. A basic rule r is of the form 

/t <— ai, . . . , a„, not bi, . . . , not bm 

and is interpreted by the function fr : 2^*°™* x 2"^*°"** 2^*°™* as follows. 

friS, C) = {h\ ai,. . . ,an e C, bi,...,bm^ S}. 

The function fr produces the result of a deductive step when applied to 
a candidate stable model S and its consequences C. 

Definition 3. A constraint rule r is of the form 

h ^ k {ai, . . . , a„, not 61, ... , not bm\ 

and is interpreted by 

fr{S,C) = {/i||{ai,...,a4nC| + |{6i,...,6^}-5| > A:}. 

The constraint rule can be used for testing the cardinality of a set of 
atoms. The rule hi ^ 2 {a, 6, c, d} states that h\ is true if at least 2 atoms 
in the set {a, 6, c, d} are true. The rule /t2 ^ 1 {not a, not b, not c, not d}, 
on the other hand, states that /12 is true if at most 3 atoms in the set are 
true. 
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Definition 4. A choice rule r is of the form 

{hi,. . . ,hk} ^ ai,. . . ,an, not 61, ... , not bm 

and is interpreted by 

fr{S,C) = {h\h£ {hi,...,hk}nS,ai,...,an G C, bi,...,bm S}. 

The choice rule is typically used when one wants to implement optional 
choices. The rule {a} b, not c declares that if b is true and c is false, then 
a is one or the other. 

Definition 5. Finally, a weight rule r is of the form 

{ai = Wai,...,an = Wa„, not bi = Wbj^,..., not bm = Wb^} > w, 
for , wi,^ > 0, and is interpreted by 

fr{S, C) = {h\ ^Wa^ + Y^ > w}. 
ai&C bi^S 

The weight rule is a generalization of the constraint rule. If every literal 
in the body of a weight rule has weight 1, then the rule behaves precisely as 
a constraint rule. 

Definition 6. Let P be a set of rules. As before we say that fp : 2"^*°™* — > 
2Atoms ig a closure if 

fpiS)= [jfr{S,fp{S)), 
reP 

and we define 

gp{S) = f]{fp{S) I fp : 2^*°™^ ^ 2^*°™^ is a closure}. 
Then, 5 is a stable model of the program P if and only if 

S = gp{S). 

The motivation for defining constraint, choice, and weight rules is that 
they can be easily and efficiently implemented and that they are quite ex- 
pressive. For example, the constraint rule 

h ^ k {oi, . . . , a„, not 61, ... , not bm} 
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replaces the program 



{/t ^ , . . . , Qii^^ , not bj^,..., not bj^^ \ ki + k2 = k, 

I <ii < ■■■ < iki <n, I < h < ■ ■ ■ < jk2 < m}, 

which contains i^'^) rules. 

Thus, a constraint rule guarantees that if the sum of the number of atoms 
in its body that are in a stable model and the number of not-atoms in its 
body that axe not is at least A;, then the head is in the model. Similarly, if 
the body of a choice rule agrees with a stable model, then the rule motivates 
the inclusion of any number of atoms from its head. A weight rule 

/l {Oi = lUaj, ... ,an = Wa„, not bi = Wb^, . . . , not bm = Wf,^} > w, 

in turn, will force the head to be a member of a stable model S if 

1^ W^ai + Wbi > w. 

Example 1. The stable models of the program 

{ai, . . . ,a„} 

false <— {ai = wi, . . . ,an = Wn} > w 
true <— {ai = ■ui, . . . , a„ = Vn} > v 

containing the atom true but not the atom false correspond to the ways one 
can pack a subset of ai , . . . , a„ in a bin such that the total weight is less 
than w and the total value is at least v. The individual weights and values 
of the items are given by respectively wi, . . . , Wn and . . . , f„. 

Example 2. The satisfying assignments of the formula 

(a V 6 V ^c) A (-.a V 6 V ^d) A{^by cW d) 

correspond to the stable models of the program 

{a, b, c, d} <— 

false not a, not b, c 
false a, not b, d 
false b, not c, not d 

that do not contain false. 
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3 The Decision Procedure 



For an atom a, let not (a) = not a, and for a not-atom not a, let 

not {not a) = a. 

For a set of literals A, define 

not (A) = {not (o) | a G A}. 

Let A~^ = {a G Atoms \ a G A} and let A~ = {a € Atoms \ not a G A}. 
Define Atoms{A) = A'^ U A~ , and for a program P, define Atoms{P) = 
Atoms{L), where L is the set of literals that appear in the program. 

A set of literals A is said to cover a set of atoms B if B C Atoms {A), 
and B is said to agree with A if 

A'^ C. B and A~ C Atoms - B. 

Algorithm 1 displays a decision procedure for the stable model semantics. 
The function smodels {P, A) returns true whenever there is a stable model 
of P agreeing with A, and it relies on the three functions expand {P, A), 
conflict {P, A), and lookahead{P, A). 

Let A' = expand {P, A). We assume that 

El AC A' and that 

E2 every stable model of P that agrees with A also agrees with A'. 

Moreover, we assume that the function conflict {P, A) satisfies the two con- 
ditions 

CI if A covers Atoms{P) and there is no stable model that agrees with A, 
then conflict{P, A) returns true, and 

C2 if conflict{P, A) returns true, then there is no stable model of P that 
agrees with A. 

In addition, lookahead{P, A) is expected to return literals not covered by A. 

Theorem 2. Let P be a set of rules and let A be a set of literals. Then, 
there is a stable model of P agreeing with A if and only if smodels {P, A) 
returns true. 
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Algorithm 1 A decision procedure for the stable model semantics 



function smodels (P, A) 
A' := expand{P, A) 
if conflict{P,A') then 

return false 
else if A' covers Atoms{P) then 

return true {A'^ is a stable model} 
else 

X := lookahead{P, A') 

if smodels{P, A' U {x}) then 

return true 
else 

return smodels (P, A' U {not (x)}) . 

function expand{P,A) 
repeat 

A' ■.= A 

A := Atleast{P,A) 

A:= AU {not x\x E Atoms{P) and x ^ Atmost{P, A)} 
until A = A' 
return A. 

function conflict{P, A) 

{Precondition: A= expand{P,A)} 
if A+nA- j^<D then 

return true 
else 

return false. 

function lookahead{P, A) 

B := Atoms{P) - Atoms{Ay, B:=BU not (B) 
while S ^ do 

Take any literal x E B 
A' := expand{P, A U {x}) 
if conflict {P, A') then 

return x 
else 

B :=B-A' 
return heuristic{P, A) . 
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Proof. Let nc{P,A) = Atoms [P) — Atoms (A) be the atoms not covered by 
A. We prove the claim by induction on the size of nc{P,A). 

Assume that the set nc{P,A) = 0. Then, A' covers Atoms{P) by El 
and smodels{P, A) returns true if and only if conflict(P, A') return false. By 
E2, CI, and C2, this happens precisely when there is a stable model of P 
agreeing with A. 

Assume nc(P,A) ^ 0. If conflict (P, A') returns true, then smodels{P,A) 
returns false and by E2 and C2 there is no stable model agreeing with A. 
On the other hand, if conflict {P, A') returns false and A' covers Atoms{P), 
then smodels{P, A) returns true and by E2 and CI there is a stable model 
that agrees with A. Otherwise, induction together with El and E2 show 
that smodels{P, A' U {x}) or smodels (^P, A' U {not (x)}) returns true if and 
only if there is a stable model agreeing with A. □ 

Let 5 be a stable model of P agreeing with the set of literals A. Then, 
fr{S, 5) C S for r G P, and we make the following observations. Let 

minr{A)= fr{C,C) 

A+CC 
A-nC=$ 

be the inevitable consequences of A, and let 

maXr{A)= \J fr{C,C) 

A+CC 
A-nC=0 

be the possible consequences of A. Then, 

1. for all r e P, S agrees with minr{A), 

2. if there is an atom a such that for all r G P, a maXr{A), then S 
agrees with {not a}, 

3. if the atom a € A, if there is only one r & P for which a € maXr{A), 
and if there exists a literal x such that a maXr{A U {x}), then S 
agrees with {not (x)}, and 

4. if not a e A and if there exists a literal x such that for some r e P, 
a G minr{A U {x}), then S agrees with {not (x)}. 

The four statements help us deduce additional literals that are in agreement 
with S. Define Atleast{P, A) as the smallest set of literals containing A that 
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can not be enlarged using 1-4 above, i.e., let Atleast{P,A) be the least fixed 
point of the operator 

f{B) = AUBU{ae mmr{B) | r G P} 

U {not a I a G Atoms{P) and for all r E P, a ^ m,aXr{B)} 
U \^not {x) I there exists a £ B such that a G maXr{B) 

for only one r € P and a maXr{B U {x})} 
U { not (x) I there exists not a E B and r E P such that 

a G minr{B U {x})}. 

Lemma 3. The function Atleast{P, A) is monotonic in its second argument. 

Proof. Observe that the function minr{B) is monotonic and that the func- 
tion maXr{B) is anti-monotonic. Hence, 

{a G minr{B) \ r G P}, 
{not a \ a E Atoms{P) and for all r G P, a ^ maXr(-B)}, 

and 

{noi {x) I there exists not a E B and r E P such that a E minr{B U {x})} 

are monotonic with respect to B. Assume that there exists a E B such that 
a E maXr{B) for only one r E P and a ^ maXr{B U {x}). li B C B' and 
a maXr{B'), then 

noi a G {not a \ a E Atoms{P) and for all r G P, a ^ maXr{B')} C f{B'). 
Consequently, both a, noi a G f{B') and therefore 

minr[f{B')^ = Atoms 

and 

maXr{f{B')) = 0. 

It follows that f{f{B')) = Atoms{P) U noi (^tom.s(P)) . Thus, is mono- 
tonic and has a least fixed point. Finally, notice that / has the same fixed 
points as □ 

We conclude. 

Proposition 4. // the stable model S of P agrees with A, then S agrees 
with Atleast{P,A). 
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Furthermore, we can bound the stable models from above. 
Proposition 5. For a choice rule r of the form 

{hi, . . . ,hk} ai, . . .,an, not 61, ... , not bm, 

let 

fr{S,C) = {h€ {hi,...,hk}\ai,...,an € C, bi,...,bm ^S}, 

and for any other type of rule, let f^{S,C) = fr{S,C). Let S be a stable 
model of P that agrees with A. Define Atmost{P, A) as the least fixed point 
of 

nB)=[jf^{A+,B-A-)-A-. 

r£P 

Then, S C Atmost{P,A). 

Proof. Note that C) is anti-monotonic in its first argument, i.e., S C. S' 
implies f'.(S',C) C f'.[S,C), and monotonic in its second argument. Fix a 
program P, a stable model S of P, and a set of literals A such that S agrees 
with A. Define 

f{B)=\Jfr{S,B) 
reP 

and 

f'{B)=[jfM-',B-A-)-A-. 

reP 

Let L be the least fixed point of /'. Since S agrees with A, 

fr{S, 5 n L) C fUA+, snL-A-)- A-, 

and f{S n L) C f'{S n L) C L. Hence, the least fixed point of f{S n •), 
which is equal to the least fixed point of /, is a subset of L. In other words, 
S CL. □ 

It follows that expand{P,A) satisfies the conditions El and E2. The 
function conflict {P, A) obviously fulfills C2, and the next proposition shows 
that also CI holds. 
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Proposition 6. If A = expand{P,A) covers the set Atoms{P) and A^ fl 
A~ = 0, then A'^ is a stable model of P. 

Proof. Assume that A = expand{P, A) covers Atoms{P) and that A'^CiA" = 
0. Then, A+ = Atmost{P, A). As fr{A+, B) C minr{A) C ^ for S C A+ , 

fr{A+,B) = f'M+,B-A-)-A- 

for every B C Thus, ^4+ is the least fixed point of 

f{B)=\Jfr{A+,B), 

reP 

from which we infer that A^ is a stable model of P. □ 
3.1 Looking Ahead and the Heuristic 

Besides Atleast{P, A) and Atmost{P, A), there is a third way to prune the 
search space. If the stable model S agrees with A but not with A U {x} for 
some literal x, then S agrees with A U {not {x)}. One can therefore avoid 
futile choices if one looks ahead and tests whether A U {x} gives rise to a 
conflict for some literal x. Since x' G expand{P, A U {x}) implies 

expand{P, A U {x'}) C expand{P, A U {x}) 

due to the monotonicity of Atleast{P, A) and Atmost{P, A), it is not even 
necessary to examine all literals not covered by A. That is, if we have tested 
X, then we do not have to test the literals in expand{P, A U {x}). 

When looking ahead fails to find a literal that causes a conflict, one falls 
back on a heuristic. For a literal x, let 

Ap = expand{P, A U {a;}) 

and 

An = expandi^P, A U {not (a;)}) . 

Assume that the search space is a full binary tree of height iJ, and let 
p = |Ap - A\ and n=\An- A\. Then, 

Op+n 
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is an upper bound on the size of the remaining search space. Minimizing 
this number is equal to minimizing 

on I OP 

log^^ = log(2^' + 2^)-(p + n). 

Since 

2inax(n,p) ^ c£i _j_ 2P <; ^^^(niV)-\-'^ 

is equivalent to 

max(n,p) < log(2"' + 2^) < max(n,p) + 1 

and 

- min(n,p) < log(2"' + 2^) - (p + n) < 1 - min(n,p), 

it suffices to maximize min(n,p). If two different literals have equal mini- 
mums, then one chooses the one with the greater maximum, max(n,p). 

4 Implementation Details 

The deductive closures Atleast(P^A) and Atmost{P^A) can both be im- 
plemented using two versions of a linear time algorithm of Dowling and 
Gallier [1]. The basic algorithm associates with each rule a counter that 
keeps track of how many literals in the body of a rule are not included in 
a partially computed closure. If a counter reaches zero, then the head of 
the corresponding rule is included in the closure. From the inclusion follows 
changes in other counters, and in this manner is membership in the closure 
propagated. 

We begin with basic rules of the form 

^ <— ai, . . . , a„, not bi, . . . , not bm- 

For every rule r we create a literal counter r. literal, which is used as above, 
and an inactivity counter r. inactive. If the set A is a partial closure, then 
the inactivity counter records the number of literals in the body of r that 
are in not [A). The counter r. inactive is therefore positive, and the rule r 
is inactive, if one can not now nor later use r to deduce its head. For every 
atom a we create a head counter a.head that holds the number of active 
rules with head a. 

Recall that a literal can be brought into Atleast{P, A) in four different 
ways. We handle the four cases with the help of the three counters. 
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1. If r. literal reaches zero, then the head of r is added to the closure. 

2. If a. head reaches zero, then not a is added to the closure. 

3. If a. head is equal to one and a is in the closure, then every literal in 
the body of the only active rule with head a is added to the closure. 

4. Finally, if a is the head of r, if not a is in the closure, and if r. literal = 1 
and r. inactive = 0, then there is precisely one literal x in the body of 
r that is not in the closure, and not (x) is added to the closure. 

Constraint rules and choice rules are easily incorporated into the same 
framework. Specifically, one does neither use the first nor the fourth case 
together with choice rules, and one docs not compare the literal and inactiv- 
ity counters of a constraint rule h ^ k {ai, . . . , a„, not bi, . . . , not bm} with 
zero but with m + n — k. A weight rule 

h<- {ai = Wai,...,an = Wa„, not bi = Wbj^,..., not bm = wj,^} > w, 

is managed using the upper and lower bound of the sum of the weights in 
its body. Given a set of literals A, the lower bound is 

aieA+ heA- 

and the upper bound is 

a^^A- b,<^A+ 

If the upper bound is less than w, then the rule is inactive, and if the lower 
bound is at least w, then the head is in the closure. 

Notice that the implementation provides for incremental updates to the 
closure Atleast{P,A) as A changes. This is crucial for achieving a high 
performance. 

Since the function Atmost{P, A) is anti-monotonic, it will shrink as A 
grows. It is no good computing Atmost{P, A) anew each time A is modified. 
Instead all atoms that might not be in the newer and smaller closure are 
found using a variant of the basic algorithm. By inspecting these atoms it 
is possible to decide which ones must be in the closure, and then the basic 
algorithm can again be used to compute the final closure. A small example 
will make the method clear. 



13 



Example 3. Suppose P is the program 

a <— b a ^ not c 

b a a <— not d, 

and suppose A has changed from the empty set to {d}. Then, we have 
already computed Atmost{P, 0) = {a, b}, and we want to find Atmost{P, A). 
If r is the rule a <— not d, then the counter of r is at first zero and then 
changes to one as d becomes a member of A. Therefore, we deduce that a 
is possibly not a part of the new closure. The basic algorithm proceeds to 
increment the counters of 6 <— a, removing b, and a b, where it stops. At 
this point the counter of the rule a not c is still zero, and we note that 
a must be part of the closure. Including a causes the counter of 6 <— a to 
decrease to zero. Consequently, b is added to the closure and the counter 
of a ^ 6 is decremented. Since nothing more remains to be done, the final 
closure is {a, 6}. 

One can argue, in this particular example, that a follows from the rule 
a <— not c and need not be removed in the first stage of the procedure. 
However, in general it is not possible to decide whether an atom is in the 
final closure by inspecting the rules of which it is a head. Notwithstanding, 
we can make improvements based upon this observation. 

For every atom a, create a source pointer whose mission is to point to the 
first rule that causes a to be included in the closure. During the portion of 
the computation when atoms are removed from the closure, we only remove 
atoms which are to be removed due to a rule in a source pointer. For if the 
rule in a source pointer does not justify the removal of an atom, then the 
atom is reentered into the closure in the second phase of the computation. 
In practice, this simple trick yields a substantial speedup of the computation 
of Atmost{P,A). 

5 Experiments 

We will search for sets of binary words of length n such that the Hamming 
distance between any two words is at least d. The size of the largest of 
these sets is denoted by A{n,d). For example, A(5, 3) = 4 and any 5- 
bit onc-error-correcting code contains at most 4 words. One such code is 
{00000,00111,11001,11110} = {0,7,25,30}. Finding codes becomes very 
quickly very hard. For instance, it was only recently proved that A{1Q, 3) = 
72 [10]. 
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Construct a program that includes a rule 

Wi not Wj^ , ■ ■ ■ , not Wj^ 

for every word z = 0, . . . , 2" such that ji , . . . , jk are the words whose distance 
to i is positive and less than d. Then, the stable models of the program are 
the maximal codes with Hamming distance d. Add the rule 

true <— m {wq, . . . , u;2" } 

and every model containing true is a code of size at least m. For the purpose 
of making the problem a bit more tractable, we only consider codes that 
include the zero word. 

The test results are tabulated below. The minimum, maximum, and 
average times are given in seconds and are calculated from ten runs on 
randomly shuffled instances of the program. All tests where run under 
Linux 2.2.6 on a 233MHz Pentium II with 128MB of memory. 



Problem 


Mill 


Max 


Average 




Problem 


Min 


Max 


Average 


^(5,3) > 4 


0.01 


0.02 


0.02 




^(6,5) > 2 


0.02 


0.03 


0.03 


yl(5,3) < 5 


0.00 


0.02 


0.02 




^(6,5) < 3 


0.02 


0.03 


0.02 


^(6,3) > 8 


0.02 


0.04 


0.03 




A(7,5) > 2 


0.05 


0.07 


0.06 


^(6,3) < 9 


0.16 


0.18 


0.17 




yl(7,5) < 3 


0.04 


0.07 


0.06 


^(7,3) > 16 


0.14 


14.19 


6.77 




^(8,5) > 4 


0.29 


0.36 


0.34 


^(7,3) < 17 


69.08 


72.29 


70.55 




^(8,5) < 5 


2.64 


2.75 


2.71 


^(8,3) > 20 


6.39 


202.41 


55.98 




^(9,5) > 6 


3.18 


8.71 


4.81 


A(8,3) < 21 


> 1 week 




A(9,5) < 7 


1127.03 


1162.10 


1145.85 



6 Conclusion 

Wc have presented some new and more expressive prepositional rules for 
the stable model semantics. A decision procedure, which has been used as 
a base for an efficient implementation, has also been described. We note 
that the decision problem for the extended semantics is iVP-complete, as a 
proposed stable model can be tested in polynomial time. Accordingly, the 
exponential worst case time-complexity of the decision procedure comes as 
no surprise. 

The literals that smodels{P, A) can branch on are, in this paper, the 
literals that do not cover Atoms{P) — Atoms (A). In previous work, for 
instance in Nicmela and Simons [6, 8], the eligible literals have also been 
required to appear in the form of not-atoms in the program. This additional 
restriction can reduce the search space, and a similar requirement is, of 
course, also possible here. The question of which literals one necessarily 
must consider as branch points is left to future research. 
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Appendix 

Let X be a set and let / : 2"^^ — > 2"^^ be a function, li A C. B implies 
/(^) C f{B), then / is monotonic. 

Lemma. Let f : 2^ ^ 2^ be a monotonic function, and let A C X . If 
f{A) C A, then Ifp (/) C A, where Ifp (/) denotes the least fixed point of f . 

Proof. Define 

S= f] A {fiX)CX). 

f{A)CA 

Then, f{A) C A implies S C A, which in turn implies f{S) C f{A) by the 
monotonicity of /. Hence, f{S) C A, and consequently 

fiS)= fl f{S)C fl A = S. 

f{A)CA f{A)CA 

Now, fiS) C S implies C which by the definition of S 

implies S C f{S). Thus, S = f{S). Moreover, for any fixed point A, 

f{A)CA imphes S C A, 

and hence Ifp (/) = S by definition. □ 

Similarly, A C f{A) implies A C gfp (/) for the greatest fixed point 
of /. Notice that if X is finite, then Ifp (/) = f^{%) for some n < \X\ 
since /(0) C //p (/). Furthermore, observe that if we are given k monotonic 
functions fi,. . . ,fk, then the least fixed point of 

k 

g{A) = U f^{A) 

i=l 

is the limit of any nest 

An+l =AnU fi(nMn), ^0 = and /j(„)(^n) ^ A„ ^ ^jfjiAn) C A^. 

In other words, the least fixed point of g can be computed by repeated 
applications of /i, . . . , /fc. 



16 



References 



[1] W.F. Dowling and J.H. Gallier. Linear-time algorithms for testing the 
satisfiabihty of propositional Horn formulae. Journal of Logic Program- 
ming, 3:267-284, 1984. 

[2] M. Gelfond and V. Lifschitz. The stable model semantics for logic 
programming. In Proceedings of the 5th International Conference on 
Logic Programming, pages 1070-1080, Seattle, USA, August 1988. The 
MIT Press. 

[3] K. Heljanko. Using logic programs with stable model semantics to solve 
deadlock and reachability problems for 1-safc pctri nets. In Tools and 
Algorithms for the Construction and Analysis of Systems, volume 1579 
of Lecture Notes in Computer Science, pages 240-254, Amsterdam, The 
Netherlands, March 1999. Springer- Verlag. 

[4] V.W. Marek and M. Truszczyhski. Stable models and an alternative 
logic programming paradigm. The Computing Research Repository, 
September 1998. cs.LO/9809032. 

[5] I. Niemela. Logic programs with stable model semantics as a constraint 
programming paradigm. In Proceedings of the Workshop on Compu- 
tational Aspects of Nonmonotonic Reasoning, pages 72-79. Research 
Report A52, Helsinki University of Technology, May 1998. 

[6] I. Niemela and P. Simons. Efficient implementation of the well-founded 
and stable model semantics. In Proceedings of the 1996 Joint Interna- 
tional Conference and Symposium, on Logic Programming, pages 289— 
303, Bonn, Germany, September 1996. The MIT Press. 

[7] C. Sakama and K. Inouc. An alternative approach to the semantics of 
disjunctive logic programs and deductive databases. Journal of Auto- 
mated Reasoning, 13:145-172, 1994. 

[8] P. Simons. Towards constraint satisfaction through logic programs and 
the stable model semantics. Research Report A47, Helsinki University 
of Technology, August 1997. 

[9] P. Simons. Smodels 2.10. http://www.tcs.hut.fi/pub/smodels/, 1999. 
A system for computing the stable models of logic programs. 



17 



p. Ostergard, T. Baicheva, and E. Kolev. Optimal binary one-error- 
correcting codes of length 10 have 72 codewords. IEEE Transactions 
on Information Theory, 45(4):1229-1231, May 1999. 



18 



